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We present a finite temperature (T) study of the t-J model on the two-dimensional triangular 
lattice for the negative hopping t, as relevant for the electron-doped NaiCo02 (NCO). To under- 
stand several aspects of this system, we study the T-dependent chemical potential, specific heat, 
magnetic susceptibility, and the dynamic Hall-coefficient across the entire doping range. We show 
systematically, how this simplest model for strongly correlated electrons describes a crossover as 
function of doping (x) from a Pauli-like weakly spin-correlated metal close to the band-limit (den- 
sity n = 2) to the Curie- Weiss metallic phase (1.5 < n < 1.75) with pronounced anti-ferromagnetic 
(AFM) correlations at low temperatures and Curie- Weiss type behavior in the high-temperature 
regime. Upon further reduction of the doping, a new energy scale, dominated by spin-interactions 
(J) emerges (apparent both in specific heat and susceptibility) and we identify an effective inter- 
action Jeff(x), valid across the entire doping range. This is distinct from Anderson's formula^, as 
we choose here t < Q, hence the opposite sign of the usual Nagaoka-ferromagnetic situation. This 
expression includes the subtle effect of weak kinetic AFM- - as encountered in the infinitely corre- 
lated situation [U = oo). By explicit computation of the Kubo-formulae, we address the question of 
practical relevance of the high-frequency expression for the Hall coefficient R'n '^ ■ We hope to clarify 
some open questions concerning the applicability of the t-J model to real experimental situations 
through this study. 
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I. INTRODUCTION 

Since the discovery of superconductivity in H2O- 
intcrcalated NCO crystals^ this compound has drawn sig- 
nificant attention. In addition, NCO displays a variety of 
remarkable thermoelectric effects, such as an unusually 
large and magnetic field dependent thermopower— and 
a T-linear Hall-coefficienfe2i&. While some of its proper- 
ties are similar to those seen in high- Tc materials there 
are some distinctive differences. Both the cobaltates 
and the cuprates become superconducting in the vicin- 
ity of the Mott-insulating state, however, while this is 
achieved through hole-doping in the cuprates, NCO re- 
quires electron-doping x between a; = 1/4 and x = 1/3^. 
While for both compounds the role of electron transport 
in two-dimensional planes - separated by thick insulating 
layers - appears to be crucial, the sites for the localized 
electron wavefunctions are organized in planar square lat- 
tices (Cu-ions) in the cuprates and in a triangular lattice 
(Co-atoms) in NCO. The most striking differences how- 
ever, are the sign of the hopping, which is inferred to 
be positive in the cuprates and negative in NCO, and 
its absolute value, which is believed to be an order of 
magnitude smaller in NCO^. While the experimental re- 
sults indicate a rich blend of phases - stemming perhaps 
from structural rearrangements with a; - it is first neces- 
sary to systematically probe the scope of a single model 
for strongly correlated systems - namely the t-J model 
- in the case of this electronically frustrated triangular 
lattice. 

In the half-filled band, where strong interactions lead 
to the spin-i Heisenberg-antiferromagnet (HAF), the 
square geometry facilitates anti-ferromagnetism (AFM) 
through its bipartite nature, while the tendency towards 



AFM in the cobaltate system is weaker. However, exten- 
sive numerical computations^ indicate Neel long-ranged 
order even in the triangular lattice. Away from half- 
filling, projected hopping on the square lattice system 
{t > 0) obeys the Nagaoka mechanism-'^''. This results in 
a competition between the ferromagnetic (FM) kinetic 
operator and AFM potential energy for finite values of 
J. Conversely, in the current case of the triangular lat- 
tice with a negative hopping t < 0, AFM tendencies are 
supported - rather than weakened - by the electronically 
frustrated kinetic energy operator. By further reducing 
the t-J model - taking J = - we explicitly isolate the 
effect of electronic frustration which leads to kinetic An- 
tiferromagnetism (kAj^ , an effect triggered only by the 
motion of holes in an AFM spin-background, even away 
from the special case of doping 1/ L, with L the number 
of sites. 

This paper has two main objectives: Firstl y, w e want 
to supplement the arguments made in ref. Illl where 
the Curie- Weiss phase in the doping range x > 0.5 was 
addressed. Secondly, we will discuss the finite temper- 
ature t-J model starting from the Nagaoka-regime and 
then extend into finite dopings. The structure is as fol- 
lows: In section II we describe our model and the exact 
diagonalization technique, which we utilize for our com- 
putations. Section III presents the results for the static 
properties chemical potential, specific heat and entropy 
as well as the spin susceptibility. We then investigate 
the dynamic Hall-coefficient Rh{ijj) and compare its dc- 
and high-frequency behavior. Lastly, we draw a con- 
nection between the Hall-coefficient and the diamagnetic 
susceptibility. For all properties, we study systematically 
the dependence on temperature, doping, and interaction 
strength. In Section IV we conclude and summarize our 
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results. 



Beginning at a doping of only a single hole in the in- 
finitely correlated system, we track the effect of kA into 
the intermediate doping regime. The specific heat devel- 
ops a low-temperature peak indicative of effective spin- 
interactions, generated by a dressed hole in an AFM-spin 
background. The same effect is manifest in the suscep- 
tibility results where we present an x- and J-dependent 
effective interaction, with a small but finite AFM value in 
the case of J = 0. We proceed by an explicit computation 
of the frequency-dependent Hall-coefhcicnt. The results 
show weak frequency-dependence of Rh and encourage 
the applicability of its high-frequency limit i?^ to prac- 
tical situations. While small in magnitude compared to 
the spin-susceptibility, we show through an identity in 
terms of derivatives - suggested by our numerical results 
- that the diamagnetic susceptibility alone may be closely 
connected to the temperature dependent Hall-coefficient. 



II. MODEL AND EXACT DIAGONALIZATION 

As a model for hole doped strongly interacting Mott- 
insulators, the t-J Hamiltonian can be applied to the 
experimental situation of electron-doped NCO after an 
electron- hole transformation^^, making use of the sym- 
metry of the Hubbard model in respect to the Mott- 
insulator. The replacement rules are t — s- — i, 6 ~ 
|1 — n| = X and —<le, where 6 is now hole-doping 

away from half-filling and n is the electron density per 
site. This leads to the t-J model 



H 



<i,j>,(7 



<«J> 



Si • Sj 



(1) 

where cj^. (cia) creates (annihilates) an electron of spin cr. 
Si is the three-component spin-operator, fii is the num- 
ber operator and i specifies the lattice site. Pq denotes 
the Gutzwiller projector and the summation is over all 
nearest neighbor pairs < i,j >. 

In the following, whenever doping x is mentioned this 
refers to electron doping in the NCO system and hole- 
doping in the model Hamiltonian. To compute the com- 
plete spectrum - as required for thermodynamic proper- 
ties - we employ the canonical ensemble exact diagonal- 
ization technique for small clusters. The computational 
demand limits the analysis to quite small systems and fi- 
nite size effects are challenging to overcome. Nonetheless, 
careful comparison between several systems of different 
geometry enables us to extract the stable behavior and 
somewhat isolate finite size effects. Conversely, finite size 
effects can help identify the favored eigenstate of a sys- 
tem. 

We have employed several toroidal and ladder geome- 
tries with L = 9, 12, and 18 sites, as well as the highly 



symmetric icosahedral geometry. The icosahedron - a 12- 
sited Platonic solid with 20 equilateral triangular faces - 
shares the property of geometric frustration with the in- 
finite system while the effect of the reduced coordination 
number (z=5) appears to have only a small, quantitative 
effect on the computational results (Fig. H^c)). Our 12- 
site torus clusters are depicted in Fig. na)-(b). A third 
12-site cluster is the simple 4 x 3-cluster (called Torus c, 
not pictured). 

While the numerical effort of exact computations for 
translationally invariant systems can be significantly 
reduced by exploiting space group symmetries of the 
Hamiltonian, these symmetries are broken upon the in- 
troduction of a magnetic field, as relevant for the evalua- 
tion of the Hall coefhcient. A magnetic field is introduced 
by the usual Peierls substitution which modifies the hop- 
ping t between sites i and j by 



t Uj (A) = t exp 



27r 



A • ds 



(2) 



where A is the magnetic vector potential and 00 the flux 
quantum. We define the flux threading a triangular pla- 
quettc as a = 1^ A • ds. In finite systems, the value of 
the smallest non-zero magnetic field is limited to values 
of a > tt/1 where I is the length of a periodic loop in 
the system. Through a particular gauge we can achieve 
I equal to the number of triangular faces in the cluster, 
this guarantees the equality of the flux values through all 
plaquettes. An example is given in Fig. Q^a). A similar 
strategy has been followed in the case of the square lattice 
quantum Hall effect We complement the computations 
on periodic systems with ladder systems, which enable an 
infinitesimal flux to be chosen due to open boundary con- 
ditions in one of the spatial dimensions (Fig. H^b)) and 
the Landau-gauge is adequate. In the case of the icosa- 
hedron (Fig. Wic)) the magnetic field is applied through 
a (fictitious) magnetic monopole in its center, creating 
equal magnetic fluxes through all triangular plaquettes. 



III. RESULTS 

A. Chemical Potential 

The chemical potential is evaluated as a derivative of 
the free energy F with respect to the particle number N, 



dF{T) 
ON 



The derivative has to be taken as a discrete difference of 
F{N). To take this derivative numerically, we use the 
central difference formula. Hence it is suitable to define 

fi+ = F{N + 1) - F{N) and = F{N) - F{N - 1) . 

The value of /i can then be obtained from 

A* — r, 
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FIG. 1: (color on-line) Small numbers along arrows (red) denote Peierls phases in units of the flux through a single plaquette, 
phases on Icosahedron (c) not shown for clarity. Dotted parallelograms (black) specify boundaries of finite clusters ((a) and (b)). 
All sites on corners of parallelograms are equivalent and translations along bounding vectors Ti and T2 tile infinite lattice; 
shaded region in (b) forbidden for ladder. Note: different coordinate systems are used for torus and ladder. For ladder-system 
we use the Landau gauge with A oc yx, (a) represents only one choice (denoted Torus a) for 12 sites. Another (denoted Torus 
b) results from ladder (b) by connecting all corners through Ti and T2 (not shown). Torus c is simple 3x4 cluster (not 
shown) . 



In Fig. Illa)-(b) we present /u(T) at several J and dopings 
X. At low temperatures /^(T) shows degenerate Fermi be- 
havior (oc T^)i^. In the strong-correlation regime ( J ~ 
and J — 0.4|t|) the ground state (g.s.) internal energy 
has a maximum at intermediate dopings 2: « 0.5 due to 
the competition of increasing number of mobile particles 
and reduction of the available vacancies upon decrease of 
doping and, hence, the low temperature chemical poten- 
tial changes sign here. In the high-temperature regime 
the doping-dependence of the Hilbert-space dimension 
determines the limit of 

^i[T -> 00) cx T 

with g{x) = (1 — x)log(2a;/(l — x)) + log(l/a;) essen- 
tially the natural logarithm of the dimension of the 
Hilbert-space. In the projected t-J model, g{x) peaks 
at X = 1/3. For a; < 1/3 the slope at large T be- 
comes negative. In Fig. |21 the comparison of /i(T) for 
different interactions shows that the impact of projected 
fermions near half-filling is to force a negative slope of 
fi for X > which occurs in the non-interacting system 
only for a; < 0. This is an example of the qualitative dif- 
ference between a strongly-interacting Fermi-system and 



its non-interacting counterpart. 



At larger x > 0.7 the effect of strong interactions is 
gradually lost and the temperature-dependence resem- 
bles that of the non-interacting Fermi-system (Fig. |2Ic)). 
For low temperatures well below the linear range of /i, fi- 
nite J > generally increases the value of the chemical 
potential with a stronger increase for x near the Mott- 
insulator. The high-temperature behavior is independent 
of J. 

In Fig. [2Jd) we compare the results with direct ex- 
perimental measurements of the T-dependent chemical 
potential"'^^ at a; « 0.83. Since our results in this doping 
range are rather insensitive to the value of U, we com- 
pare with our U — results. Here, IJ.{T) was obtained 
directly from the tight-binding density of states by the 
implicit equation (n) —2x f(e, fi)p(e)de with fixed 
(n) = 1 + X, where / denotes the Fermi-Dirac distri- 
bution function and p is the non-interacting density of 
states. This fit confirms the notion of NCO as a system 
with an unusually small Fermi temperature scale. In Fig. 
|2Id) we show as a shaded region the values t = —100 
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FIG. 2: (color on-line) Chemical potential /i(T) for all accessible dopings on Torus h: x — 2/12 = 0.16 to 10/12 = 0.83 for J = 
(a), J = 0.4\t\ (b) and C/ = (c). Data for [7 = computed for infinite non-interacting system, comparison with (a)-(c) shows 
relatively weak effect of interactions in the doping range x > 0.8.(d) Squares are experimental data for Afi{T) — ^jl{T) — /i(0) 
from ref. [Tsl Shaded area is a comparison with our data for (7 = for t = —100 K to t = —300 K. 



K to t = —300 K, photoemission results^ point to a 
value of \t\ K, IQmeV ^ —118 K, hence |i| here is more 
than an order of magnitude smaller than in the cuprate- 
superconductors. In the following, when making contact 
with experiments, we take t = —100 K for simplicity. 



B. Specific Heat, Entropy 
In Fig. ISta)-(b) we present the specific heat per site 



1 d{E) 
L dT 



with {E) the thermodynamic expectation value of the in- 
ternal energy. We first concentrate on the case of J = 0, 
corresponding to an infinite on-site interaction. Near the 
band insulator (a; « 1) the temperature-dependence re- 
sembles that of non-interacting fermions, computed at 
fixed particle density by first determining fJ,{T) and then 
the T-derivative of the internal energy. When the doping 
is reduced, the main energy scale - manifest in a peak of 
Cv{T) at temperature Tm{x, J) - continuously moves to 
smaller temperatures, away from the tight-binding peak. 
Even at our lowest doping x = \/L the peak-feature per- 
sists (called kA-peak). This feature is magnified in Fig. 
21 for several different systems. 



Numerical work on the g.s. of the problem of a single 
hole in an AFM spin-backgrounciS, suggests this system 
to exhibit a weak effective AFM interaction J^/f (x 1/L 
due to the projected kinetic energy term only. While 
for one hole this excitation decays in the thermodynamic 



limit, J^ff nonetheless induces a Neel ordered AFM g.s. 
and low energy excitations similar to those observed in 
the HAF^. 

The validity of our computations for the case of a single 
hole is underscored by the entropy. In Fig. the tem- 
perature dependence of the entropy shows that the low- 
energy features observed in the specific heat are necessary 
for the entropy to obtain its (known) high-temperature 
limit. For T/\t\ w 0.075 half the entropy is recovered 
{s{T)/s(oo) = 0.5), indicative of entropy-enhancement 
due to electronic frustrationiSi. For intermediate dopings, 
the entropy reaches half its maximum at T/\t\ 0.75 and 
for dopings x « 0.7 - relevant to the Curie- Weiss phase 
of NCO - we find T/\t\ « 0.78. 



Within a high-temperature expansion^^ a numerical 
argument was given for an S-shaped low temperature be- 
havior of the entropy for the case of x = 0.2 and x — 0.4. 
Our computations at small dopings rather suggest s ~ T 
for low T as applicable to a metal from Fermi-liquid ar- 
guments. It is remarkable here that Cy{T) refiects these 
low energy excitations, notably even for larger dopings 
than the one applicable to the Nagaoka-situation. 

We now distinguish the effect of J > (Fig. Ofb)). 
This scenario has been investigated by several authors 
for the case of the hole doped square lattice and the 
HAFiiLi^. For the situation of a single hole J > shifts 
the fcj4-peak to larger temperatures Tm — J/t. In Fig. 
IHtb) we present the peak positions as function of doping 
and J. We find Tm{x, J) to increase monotonically both 
with X and J and the slope as function of J for fixed x 
to decrease for larger values of x. Extrapolating these 
results at a given value of J to x = we obtain roughly 
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(a) T/|ti (b) T/|ti 



FIG. 3: (color on-line) Specific heat per site as function of T for varying x (computed on 18 site {x > 0.72) and 12 site cluster 
{x < 0.72), (a) J = 0, (b) J = 0.4|t|. Dotted lines show Cv{T) for non-interacting fermions and solid line for exact results. 
Note: scale changes in plots of different dopings, long temperature exponential behavior at large x due to finite-system induced 
gap in energy spectrum. Results for x — 0.72 appeared in ref. ^3 




' 1/18, Torus 
1/14, Torus 
1/12, Torus a 
1/12, Torus b 
1/12, Torus c 
1/13, Torus 
1/9, Torus 
1/12, Icos. 



0.2 0.3 

T/|t| 



0.4 



0.5 



FIG. 4: (color on-line) Cv{T) at x = 1/L for several sets 
of boundary conditions, peak in specific heat is system- 
independent property. 



Tm{x = 0) ~ J. Our result should be distinguished from 
investigations on the t > Q square lattice^^, where T,„ vifas 
found to follow the opposite trend decreasing slightly with 
X at small dopings, along with a rapid suppression of the 
peak. This supports the notion of electronic frustration 
contributing to AFM correlations in the present study 
in opposition to competing AFM exchange coupling and 
the FM Nagaoka-mechanism in the case of t > 0, capable 
of quickly eliminating the low-energy spin excitations. 



An experimentally interesting quantity related to the 
specific heat is 7 defined by Cy{T) ~ 7T or 7 = ds/dT 
corresponding to the low energy electronic part of the 
specific heat. When extracting this quantity from our 
numerical results, we neglect the contributions at very 
low T arising purely from the finite system induced 
gap in the spectrum (scaling as exp(— A/T) where A is 
the gap between the g.s. and the first excited state). 
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FIG. 5: (color on-line) Entropy per site for three different 
dopings at J = (x = 1/12 = 0.08, x = 5/12 = 0.42, 
a; = 13/18 = 0.72), horizontal lines indicate high-temperature 
limit s{T — > oo) = kB\og{g{x))/N, x = 0.08 corresponds to 
Fig. 0] Exponential increase at low T for intermediate x due 
to finite system-size gap A. 



For the large doping case x — 0.7 we foundii 7 « 
45 (±5) mJ/molK^ (70 « 33 mJ/molK^ for non- 
interacting tight-binding electrons), indicating a signif- 
icant many-body enhancement. Experimental resultsiS 
indicate 7 « 26 mJ/molK^. On the other hand, in the 
low-doping regime the linear increase with T is only ob- 
served in a very small temperature range and this range 
is possibly better classified as a non-Fermi-liquid regime. 
Roughly, 7(x)/7o(a;) ^ 1 for a; > 0. The increase in 
7 with decreased doping is also demonstrated clearly by 
the low T slope of s(T) in Fig. 



C. Spin Susceptibility 

To gain further insight into the role of spin-fluctuations 
and to address the existence of the experimentally ob- 
served Curie- Weiss phase of NCO^S^ we compute the spin 
susceptibility xi'^) function of temperature from the 
Gibbs free energy G{T, M), here M = J2^ Sf is the total 
magnetization. This quantity follows from the Helmholtz 
free energy F(T, H) through a Legendre-transformation 
G{T, M) = F{T, M) - H§^ = F{T, H) + MH. H is the 



externally applied magnetic field. 



X 



can then be evaluated numerically, by constraining the 
canonical ensemble to the sector of minimal 15*°*!, and 
those sectors where 5'*°* is raised (lowered) by h. 



For low temperatures, near the band-limit [x = 1) 
we identify weakly temperature-dependent, Pauli-like be- 
havior of X, reminiscent of the experimentally observed 
spin density wave (SDW)-phase (Fig. O. When x is low- 
ered past X — 0.75 we find^^ a crossover to a Curie- Weiss 
like behavior with appreciable renormalizations of the g.s. 
susceptibility compared to its bare value, mirroring the 
experimental findings. This behavior can roughly be un- 
derstood by picturing electrons distributed on an empty 
lattice: Initially, the many-body eigenstate enables the 
electrons to move around without interacting directly 
with their neighbors, yielding a susceptibility that is J- 
dependent only to second-order. When the density is in- 
creased {n > 0.25 or equivalently n < 1.75), electrons are 
forced to interact through the exchange-coupling directly 
and hence yield a strongly J-dependent magnetic suscep- 
tibility. The emergence of the Curie- Weiss phase is fur- 
ther supported by the strong magnetic field suppression 
of the Heikes-Mott term in the thermopower within the 
t-J model, an effect due to spin-entropy^ as was demon- 
strated in ref. 

For temperatures larger than a certain crossover tem- 
perature (which we call T*), for all dopings, x{T) can be 
described through the Curie- Weiss form 



xiT) 



1 L 



with the Weiss-temperature and the effective magnetic 
moment Pef / , is the volume of the crystal and L the 
number of sites. For all accessible system sizes - in par- 
ticular Torus a,b,c (Fig. ^) - we compute x(T) to ex- 
tract 9, Peff and T*, the parameters characterizing the 
magnetic properties at a given doping and interaction 
strength. To determine 9 and Peff we fit the high- 
temperature part of our computations for x(r) to the 
function f{T)^a/{T-b) (compare Fig. Cfa)). In Fig. 
[7fa)-(c) we present the inverse susceptibility at several 
different dopings near the Mott-insulating phase. The 
results for 9 and Pe// are shown in Fig. [7fd)-(f). For 
all X and J > 0, 9 remains negative indicating AFM cor- 
relations. 9 shows a linear increase when the doping is 
increased or the interaction is decreased. In short, as 
mentioned in ref. IllL 9 can be collapsed into the form 
9{x,J) — —c Jeff{x) where 

Jeffix) = J(l -I- c"x\t\) + c'x\t\ 
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FIG. 6: (color on-line) Spin susceptibility per site xiT) for all accessible dopings x on 12-site torus. Dependence on J/\t\ —{0, 
0.1, 0.2, 0.3, 0.4, 0.5), indicated by arrow, blue dotted line shows bare-susceptibility for non-interacting fermions of density 
n = 1 + x and black dotted line indicates susceptibility of isolated spin-i fermions of density n = 1. Plots for x > 0.58 appeared 
in ref. 11. 



with c = 4.0, c' = 0.01425 and c" = -0.9175. In contrast 
to the relation given by Anderson et alA in the case of 
the t > square lattice, in our situation of i < even 
at J = there is a slightly negative intercept suggesting 
AFM spin-interactions in the low temperature regime of 
the electronically frustrated U = oo Hubbard model, in- 
dicative of kinetic Antiferromagnetism (kA ). 

On the other hand, once J > is allowed, this inter- 
cept grows strongly and expresses dominant AFM corre- 
lations (Fig. E{e)). We compare these results with ex- 
periments for various dopings [x > 1/2) of NCO. Among 
others, the cases of x = 0.75^W^, x ^ 0.7^2, x = 0.6 
and x = 0.9^^1, various x in the three-layer system'^^ 
and X — 0.55^^ were experimentally investigated. The 
comparison with experiment proves difficult due to the 
large number of structural phasesi^ij the role played by 
the Na-ions in the crystal and the large quantitative and 
qualitative difference between powder and single crystal 
samplesS^. However, the Weiss-temperature appears to 
lie between 100 K and 200 K near x = 0.7 (varying con- 
siderably with H\\c and iJ ||a^'^) and was measured to be 
smaller at larger dopings^^. Hence, at least for x > 0.5, 



the t-J model description yields reasonable experimental 
agreement. 

The effective magnetic moment Peff iu units of /is 
has a value close to that obtained for localized spin-i 

particles at half- filling 111/2 = <? x s{s + 1) = \/3 and 
decays linearly for x ^ 1. 

The inverse susceptibility shows a minimum at a tem- 
perature T* which becomes more pronounced with in- 
creasing J. At low temperatures the system prefers an- 
tiferromagnetic order which is strengthened by J > 0. 
Thermal spin fluctuations break this order and when T is 
increased to values of order J a small perturbative B- field 
is sufficient to align spins antiparallel to the field. Hence, 
the susceptibility is largest at this temperature. When T 
is further increased, the ability of the field to flip spins 
is weakened by thermal fluctuations, consequently, x(T) 
decreases. The behavior of this peak as function of hole 
doping away from the Mott-insulator has already been 
investigated within a high-temperature expansion^^. 

Due to the small size of the systems available to our 
study, it is difficult to confirm the results for the doping 
dependence^^. Here, we study the behavior as function 
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FIG. 7: (color on-line) (a)-(c) Inverse susceptibility vs. temperature for dopings x — 1/12 = 0.08 (bottom), x — 2/12 — 0.16 
(center), and x = 3/12 = 0.25 (top) for Torus a, b and c (1 to r). Solid vertical arrows show increase in J, dotted (blue) 
line shows example for fit to extract Weiss-temperature 9{x,J) and effective moment Peff- Dotted vertical arrows indicate 
temperature T*; Note: comparison of the three systems allows isolation of finite size effects, most visible for x — 0.16. (d) 
Effective magnetic moment Peff- (e) Weiss-temperature 9 vs doping x for J — . . . 0.5. Letters in angled brackets correspond 
to experimental data from [a] Viciu— et al., [b] Gavilanc'^^ et al., [c] Sales;^^' et al., [d] Takeuchi-A et al., [e] Motohashi^^ et al., 
and [f] Effective moment of isolated spin-|-particles on half-filled lattice, (f) Magnification of 6 vs. x for J — for several 
clusters. 



of J, considering several dopings near half-filling and sev- 
eral system geometries (Fig. |HIa)). The figure indicates 
that T*/|i| roughly scales as J/\t\, hence, the scale at 
which the system undergoes a transition from an AFM- 
ordered to a paramagnetic state is set by the interaction 
strength. As was shown for the case of J = and t > 
in the hole doped system^, a single hole favors AFM Neel 
order and numerical results^ indicate that the triangular 
lattice Heisenberg model also favors Neel order - at least 
in the g.s. and the elementary excitations - hence it is not 
surprising that a finite J > emphasizes this tendency. 
This is in opposition to the case of < < where the ki- 
netic energy term favors ferromagnetism and competes 
with the Heisenberg-interaction preferring AFM. 



D. Hall Coefficient 

Experimentally, NCO has been shown to exhibit un- 
usual thermoelectric properties, resulting from strong 
electron correlations. Wang et al.^ measured the dc-Hall 
coefficient as function of temperature at several dopings 
showing unexpected temperature dependence with an os- 
cillatory behavior below 200 K and a rapid linear increase 
for temperatures above 200 K. 

Much theoretical work has been in- 
vested into the study of Rh on square 

^^^j-^„„ j3. 12. 27. 28. 29. 30. 31. 32. 33. 34. 35. 36. 37. 38. 39 both for 

g.s. properties and for finite temperatures. For 
the triangular lattice, it was possible to predict the 
high-temperature behavior within a high-temperature 
expansion of the t- J-model^^^^ in its high-frequency 
limit. Several authors have also attempted to describe 
this behavior through Boltzmann transport theory^Si^ 
within the relaxation time approximation, where Rh is 
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FIG. 8: (color on-line) (a) AFM-paramagnetic crossover- 
temperature as function of J for several dopings on 12-site 
torus clusters a, b, c. (b) Evolution of T,„ as function of J, 
values averaged over three 12-site clusters a, b, c. 



expressed as an integral over the Fermi-surface curvature 



and the relaxation time is taken to be isotropic. 
For large dopings - x > 0.7 - where the hole- like 
Fermi-surface of NCO approximates a circle, this leads 
to the same high-temperature linear coefficient as the 
t- J-high-temperature expansion, however, the behavior 
completely diverges in the dense regime (x < 0.7), 
possibly due to the misleading replacement of the inter- 
acting Fermi-surface by its non- interacting counterpart. 
Also, at low temperatures - where the relaxation time is 
anisotropic - the oscillatory behavior seen in experiments 
is not reproduced. 

To shed some light on the low-temperature and 
strongly-correlated behavior, we start from the high- 
frequency limiti for Rh{T,uj): 



R*h{T) = lim Rh{T,uj) 

LO — >OQ 

ivL < [Jx, Jy] > 



lim 



(3) 



Here, e is the electron charge, Jx {Jy) are the currents 
in X (?/)-directions, v is the unit cell volume, and t^.x 
is the diagonal part of the stress tensor. To ensure the 
effect of strong correlations within this limit, it is neces- 
sary to force the hierarchy t < uj <^ U . We enable this 
by projecting out double occupancy using the Gutzwiller 
projector, hence the t-J model is the appropriate descrip- 
tion. 

Our plan is as follows: we show the connection of this 
high-frequency limit - not measurable in experiments - 
to the dc or low-frequency behavior. Then we evaluate 
the dependence on the parameters T, J and x for several 
systems. Hence, for the first part we compare i?^ with 
the expression derived from the Kubo formulae for the 
electrical conductivities 



^ 1 exp-(/3ei.) - exp-(/3e^) \ 
ujv y Z — — U! — irj J 



with ei, the v'th energy eigenvalue of H and Z the canon- 
ical partition function. In terms of this expression the 
Hall coefficient is given by 

rhh=v^( ■ (5) 

oB \axA^)^y,y{^) J b=o 

As a transport object Rh is a demanding quantity to 
study on finite clusters. The effect of boundary con- 
ditions and their finite size implications are difficult to 
eliminate. We introduce a small magnetic field through 
the Peierls substitution (Eq. |2Jl which impacts only on 
the orbital motion of the electrons and does not couple 



to the spin. Studying Rh on toroidal and spherical ge- 
ometries reduces the finite size effects on the transport, 
however, these periodic systems constrain the smallest 
accessible flux 0o to tt/L. For L = 12 this is already a 
considerable flux. An alternative are ladder systems such 
as in Fig. ^b) where the limit (/)o ^ is possible. How- 
ever, ladders have additional finite size effects impacting 
on the transport quantities. 

For several dopings (x = 0.08, 0.67, 0.75, 0.83) we com- 
pute Rh{uj) as a function of temperature both on the 
torus and the icosahedron. To take into account the fi- 
nite size of the system, a broadening i] of the energy equal 
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to the mean level spacing of the current matrix elements 
was employed, this requires 77 « 3|t|. In Fig. |^we show 
Rh as function of temperature at three different frequen- 
cies. Comparing a; = to w = 00 we find that the de- 
limit is qualitatively the same as the high-frequency limit, 
however, the overall temperature-dependence is scaled by 
a renormalization factor s = Rh{TtUJ = 0) / Rh{T, 00) 
where s is only weakly temperature dependent. For 
the torus, s was found to be 0.83 (0.78) at a; = 0.83 
[x = 0.75). At X = 0.083 the many-body renormaliza- 
tion is stronger. 

Conversely, in Fig. ^| we show the frequency- 
dependence of Rh at several temperatures and different 
dopings. Similar to the work done by Grayson et al4£ 
on the Hall angle, finite frequency experiments may be 
helpful to confirm the weak frequency-dependence seen 
in our computations: With \t\ w IQmeV results in the 
frequency range lo > 3|t| ~ iOmeV could be compared 
to our theoretical results and hence lend experimental 
support for the usefulness of the high-frequency limit^ of 
Rh- 

We have computed two additional quantities related 
to the dc-limit as checks of the validity of the dc- 
computations. Zotos et al*^ derived an expression for the 
g.s. expectation value, Rh{T = 0) = tjt ^ log D{n) with 



D the charge stiffness D 
and 



the g.s. energy 

a constant flux inducing a persistent current^. In 



this fashion, we evaluated the doping dependence of the 
g.s. and found R^ > in the experimental regime near 
x « 0.7 and a divergence i?^ ^ — 00 as the Mott-limit 
is approachedii. The second check refers back to the ba- 
sic definition of the dc-Hall-coefficient Rh = —Ey/{Bjx) 
and computing this directly on a ladder system. This 
computation showed no change of sign for x ~ 0.7 as 
function of T and a qualitatively similar behavior to the 
Lu — 00 result. 
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FIG. 9: (color on-line) Rh{ui,T) on 12-site torus at several 
frequencies for several dopings and J = 0. 
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FIG. 10; (color on-line) Rh{uj) at several temperatures and 
dopings as function of frequency on 12-site Torus a for J — 0. 
Temperatures in all plots T/\t\ = (0.5, 1, 2, 4). 



Encouraged by the notion of weak frequency- 
dependence, we studied i?^ = Rh{^ = oo). We have 
studied several 12 site tori, several doping values on an 
18 site torus, the icosahedron and several ladders. 

The g.s. R*h{x) shows a Drude drop-off near the band 
insulator and a saddle point around x = 0.5 as was 
reported by the authors in ref. Ilj When the strong- 
correlation regime near the Mott-insulator is approached 
(x < 1/2), the slope again increases and for 1 hole 
(x = 0.083) R% changes sign. 

The high temperature behavior confirms the prediction 
from the high-temperature series expansionSiiS. The in- 
teresting temperature range is for T < 2\t\ where the 
experimental data shows an oscillatory behavior and an- 
alytical results are not available (Fig. Ill() . Our com- 
putations show a distinct minimum in Rh{T) which ap- 
pears near T = \t\ in the Curie- Weiss regime {x « 0.7) 
and continuously moves to smaller temperatures as the 
doping is decreased. For the smallest dopings accessible 
{x = 1/L) this minimum appears near T = 0.2|t|, while 
the behavior becomes linear at higher T in all cases. In 
the Curie- Weiss regime some systems also show a weak 
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maximum close to T = 0.25|t|. These results point to- 
wards the existence of two distinct energy scales. The 
pronounced minimum appears to be correlated with the 
peak in the electronic specific heat (Fig. Ill|) . which also 
occurs around T = \t\ in the band limit, and gradually 
shifts towards smaller temperatures with decreased dop- 
ing. 

To strengthen this point it is useful to include the de- 
pendence on J. As reported for the specific heat, the 
effect of J is to move the peak in specific heat to larger 
temperatures. This behavior is mirrored by Rh- In Fig. 
Illlwe show Rh{T) for several dopings and J. At x = 0.75 
the peak of Cy{T) shifts from T = 1.3\t\ to T = 1.5\t\ 
while the minimum in Rh{T) moves from T — 0.75|i| to 
T = l.l\t\. At X = 0.67 the behavior around T = \t\ is 
similar, however an additional peak in the specific heat 
emerges at T < 0.5|t|. A qualitatively similar feature is 
seen in other clusters, however, its magnitude is depen- 
dent on the system geometry. Hence, we do not attempt 
to extract quantitative information in the T-range below 
T = 0.5|i|. Nonetheless, it is important to note that 
with decreasing doping a second energy scale - caused by 
spin interactions - becomes dominant. At small doping 
a; = 0.16 and x = 0.08 the specific heat becomes strongly 
J-dependent and the main peak in the specific heat - 
as discussed in the first section - increases significantly 
with J. The variations in Rh(T) with J are also more 
pronounced. 

Assuming t « —100 K we compare our results to the 
experimentally available datafe at a; = 0.31 and x = 0.71, 
this was done in ref. for x = 0.71. In Fig. 1121 the 
comparison at a; = 0.31 shows good agreement for the 
slope at high T and reproduces the temperature scale for 
the minimum near 25 K. For x = 0.71 the temperature 
scale for the minimum near 100 K is seen both in ex- 
periment and theory. The slope of the high temperature 
limit differs by a factor of 4 and was adjusted in the plot 
for better comparison. Also, the theory does not show 
the change of sign as function of temperature. However, 
the oscillatory behavior below T = 100 K and the onset 
of the linear increase for T > 200 K agrees well with ex- 
periment. Hence, the exact temperature dependence of 
Rh remains a challenging task - both on the theoretical 
and experimental side. 



E. Diamagnetic Susceptibility 

Lastly, we study the orbital susceptibility per site 

where F is the free energy and B the magnitude of the 
perpendicular magnetic field. We define the dimension- 
less fiux a = BA/(f>Q through a triangular plaquette with 
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FIG. 11: (color on-line) (left) R*u{T) for several x and J for 
Torus b. (right) Cv{T) for several x and J on same Torus. 



A = V3aV4 the area of the plaquette and a the lattice 
spacing. We get 



1 d'^F 



Xd = --r 



L da^ 



da 
dB 



It is useful to employ ladder systems as shown in Fig. 
^b). Defining t = h'^/2mP we have 



1 



2to 2 2 1 9f 



L da 



where / is the dimensionless F/t and is the Bohr- 
magneton. In Fig. I13r a) we present the diamagnetic 
susceptibility and its first two derivatives with respect to 
T . For free electrons in an infinite system, Xd corresponds 
to electrons moving in circular orbits with radii inversely 
proportional to the magnetic field. Rh measures the po- 
larization induced in a system with a stationary current 
and perpendicular magnetic field. It is hence reasonable 
to ask for a connection between these two quantities. Our 
numerics suggest the relation 



T 



d^Xd 



c{x) 



(6) 



where c(x) is a function depending on doping x. 

In Fig. I13r b) we display our numerical result 
for R*fj{T) at a doping relevant to the experimental 
situation^ compared to the result obtained from integrat- 
ing Eqtn. @ and choosing the two integration constants 
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FIG. 12: (color on-line) Comparison of 12-site Torus b and 
icosahedron with experiment'^'' at x = 0.31 and x = 0.71, we 
take t — —100 K, in the lower plot our numerical results are 
scaled to match the slope of the linear part at T > 200 K in 
the experimental results. 
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appropriately. For J = this was done by the authors 
m ref. in Ignoring the range for T < 0.5 |t| which may 
be strongly influenced by finite size eS^ects, we find that 
both curves show a minimum at T « \t\ and then increase 
until T = 2\t\. The reaction to finite J > is also sim- 
ilar for both curves, moving the minimum from T = \t\ 
to roughly T = 1.25|t|. We plan to further investigate 
this connection within a high-temperature expansion for 
XdiT), in particular for the doping regime close to the 
Mott-insulator. 



IV. CONCLUSION 

In summary, we have presented a systematic study of 
the entire doping regime for the triangular lattice t-J 
model with electronic frustration (t < for n > 1). We 
have identified low energy excitations, stemming purely 
from kinetic processes (J — 0) of infinitely correlated 
spin-i fermions. These are manifest in a peak structure 
in the specific heat emerging at very low T near the Mott- 
insulator, however persisting for intermediate dopings. 
The effect of electronic frustration is further supported 
by entropy-enhancement in the underdoped regime with 
a low entropy recovery temperature. 



FIG. 13: (color on-line) (a) XdiT) and unitless Hall-coefRcient 
vh and their derivatives on a 12-site ladder at a; = 0.83. (b) 
Comparison of numerical results for x = 0.83 for R*u{T) (ob- 
tained on an 18-site torus) with the integral of Eqtn. Q 
{Xd(T) obtained on a 12-site ladder srotem) for two different 
values of J. Appeared partly in ref. till 



We have studied both the diamagnetic and spin- 
susceptibility. For the spin susceptibility we extract a 
relation for effective spin-interactions 9{x^ J) valid for all 
dopings. This stands in contrast to Anderson's formula 
for FM < > 0, since in our case AFM is supported rather 
than suppressed by doping. In the experimentally rele- 
vant doping range of a; > 0.5 we offer an explanation for 
the crossover between the Curie- Weiss and SDW metallic 
phases, as observed in experimental. For the diamagnetic 
susceptibility Xd a connection is obtained between Xd and 
the temperature dependent Hall coefficient. 

The dynamic Hall coefficient reveals a weak frequency- 
dependence and a qualitatively similar dc and high- 
frequency behavior, with stronger renormalizations close 
to the Mott-limit. This strengthens the practical 
relevance of the high-frequency limi^ R*fj. The 
temperature-dependence confirms the linear slope at high 
temperatures^ and reproduces the experimentally ob- 
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served minimum near T=100Kata; = 0.71. 

In conclusion, the t-J model is applicable to NCO for 
X > 0.5 and shows that the interface of strong corre- 
lations and metallic behavior indeed describes the ex- 
perimental situation here. For x < 0.5 the frustrated 
t-J model leads to subtle effects of strong correlations, 
even in the case of J = 0. Our results show, that for 
the proper description of experimentally observed normal 
state properties^ there - such as the charge-ordered insu- 



lator at a; = 0.5 or the paramagnetic metal for x < 0.5 - 
an extension of the simple t-J model would be necessary. 
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